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The ab-initio calculation of resonance lifetimes of metastable anions challenges modern quantum- 
chemical methods. The exact lifetime of the lowest-energy resonance is encoded into a complex 
"density" that can be obtained via complex-coordinate scaling. We illustrate this with one-electron 
examples and show how the lifetime can be extracted from the complex density in much the same 
way as the ground-state energy of bound systems is extracted from its ground-state density. 



INTRODUCTION 

Shape and Feshbach resonances have been observed in 
many molecules both in the gas-phase and on surfaces. 
These states often result in the formation of molecular 
negative ions whose energies and lifetimes can be mea- 
sured experimentally pQ. Due to their diffuse electron 
clouds, finite lifetimes, and subtle manifestations of elec- 
tron correlation effects, the calculation of the electronic 
structure of such metastable anions poses special diffi- 
culties for the theorist (see recent review by Simons [5]). 
Obtaining accurate resonance lifetimes from first princi- 
ples is a formidable task both because of the many-body 
nature of the problem (in particular bound-free correla- 
tion [3J) and because these states occur in the energy 
continuum where no normalizable wave function exists. 

Standard Density Functional Theory (DFT) gj [5] 
cannot come to the rescue. Even the exact exchange- 
correlation (XC) functional of DFT does not predict the 
true values for the anions' negative electron affinities be- 
cause the correct solution converges in the infinite basis- 
set limit to the ground state of the corresponding neutral 
species (plus one electron off to infinity) , and not to the 
metastable state of interest. Artificially binding the run- 
away electron with the use of a finite set of localized basis 
functions yields reasonable (approximate) resonance en- 
ergies [Bj. Extrapolation techniques also work well for 
this purpose [7], but the lifetimes are simply not accessi- 
ble via standard DFT. 

Time-dependent DFT (TDDFT) is needed as a 
matter of principle for high-lying resonances (e.g. see 
ref. [9] for autoionizing resonances), but a complex- 
scaled version of ground-state DFT may be enough for 
the lowest-energy one. After all, in many respects this 
resonance plays the role of the 'ground state' of the 
metastable anion [10]. 

This letter reports our proof-of-principlc demonstra- 
tion that the 'ground-state' resonance lifetime is encoded 
in the complex resonance density, and that the lifetime 
can be extracted from this complex density in much 
the same way as the ground-state energy is extracted 
in bound systems from the ground-state density. 

Complex- Scaling: For few-electron systems or simple 
model systems, the complex-scaling method [HJD21H3] ls 



well-established and has drawn much attention since the 
original theorem of Balslev and Combes [13] ■ It has been 
applied successfully to study resonance phenomena in 
many branches of science: atomic collision processes |15j . 
acoustic resonances in tunnels |16] „ quantum-confined 
Stark effect in quantum wells |17j . and many others. The 
method has been formulated to treat the many-electron 
problem within mean-field theories like Hartree-Fock [18 
as well as straightforward configuration-interaction ex- 
pansions |19j ; it has also been used in combination with 
Density Functional Theory to study the broadening of 
levels in the vicinity of metal surfaces j20j [21] [22], or 
even with TDDFT to study above-threshold ionization 
of small negative ions [23]. The fundamental question of 
scaling of density functionals upon coordinate-rotation in 
the complex plane, however, has not been investigated. 

Here are the most relevant aspects of the complex- 
scaling technique [T21 121] : When all coordinates of a 
many-electron Hamiltonian H are multiplied by a phase 
factor e , the resulting non-hermitian operator Hg has 
the following properties for a suitable range of 6 : (1) a 
discrete spectrum of real eigenvalues that is identical to 
that of the original (unrotated) Hamiltonian; (2) Scatter- 
ing thresholds that are also identical to those of the orig- 
inal Hamiltonian; (3) continua beginning at each scat- 
tering threshold, rotated by an angle 26 with respect to 
those of the original Hamiltonian; and (4) Complex dis- 
crete eigenvalues appearing in the lower-half of the com- 
plex energy plane. Their corresponding eigenfuctions are 
square-integrable and can be associated with resonances 
whose energies and lifetimes are determined respectively 
by the real and imaginary parts of the eigenvalues. 

COMPLEX DENSITIES 

We focus attention on the complex density ng(r) asso- 
ciated with the lowest-energy resonance, 

ng(r) = {^\h(r)\^f) (1) 

where n(r) is the density operator, and (^g\ and l^fg') 
are the left and right eigenvectors of Hg corresponding 
to the lowest-energy resonance. We require that ng(r) be 
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normalized to the number of electrons, as real densities 
are: 



J drn e (r) = N , 



(2) 



and ask whether ng(r) contains the information about the 
lifetime of the resonance. Our hypothesis, prompted by 
partial empirical evidence [S] and previous observations 
PHI [25] , is that the energy and lifetime of this resonance 
are determined respectively by the real and imaginary 
parts of the ground-state energy functional E[n] evalu- 
ated at ng(r). More specifically, the complex E[ng] is 
equal to: 



£[n e ] - -C-^ne] = F g [n e ] + J drn e (r)v(ve i6 ) (3) 

where £ is the lowest resonance energy, C the corre- 
sponding lifetime, v the external potential, and Fg the 
complex-scaled Hohcnbcrg-Kohn universal component of 
the energy density functional jl]. 

Our purpose here is to test the validity of Eq.(j3]) for 
one-electron systems, where Fg[n] is known exactly; it is 
given by F g [n] — e~ 2%e T w [n], where T w [n] is the Von 
Weizsacker functional |2"6"1. 



Tw[n] = g J dr 



lV^(r)| 2 
n(r) 



(4) 



Analytical Test: First, consider a radially-symmetric 
external potential given by v(r) = — \r~ 2 + r _1 . This 
is admittedly not too chemically meaningful but has the 
nice features that its resonances are known analytically 
[2"T] . It has no bound states but a resonance of energy 
£ = 1/4 and lifetime L — 2/y/3. Its complex-scaled 
Hamiltonian 
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(5) 



has the following radial eigenstate: 
4> e {r) = N e r^- i ^' 2 e irk6 ; kg = -ie l6 (l + iV3)/2 (6) 

The function 4>g{r) is not square-integrable for arbi- 
trary 9, but in the range [tt/6 < 9 < 7ir/6] it is, and 
the constant Ng can be chosen in this range as Ng — 

- 1 1/2 

(— 2ikg) 2 ~ l ^3 /Y{2 — iy3) , so the normalization con- 
dition for the resonance density, Eq.(j2|, is satisfied. Since 
the adjoint of the hamiltonian (|5| is simply equal to its 
complex-conjugate, the radial resonance density is given 
by hg(r) — j\fZ r 1 ~' l ^e 2 ' lker (we use a tilde on top of 
n to denote it is the radial density, h(r) = 47rr 2 n(r)). 
Even though this function lies outside the domain of real 
densities for which the functional F[n] was established 
[3], we plug it nevertheless into the right-hand-side of 



Eq.(|) with Fg[n g ] = e~ 2ie T w [ne], and find the result: 
E[ng] = j — i^-, whose real and (minus 2 times the) 
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FIG. 1. The bound electron density of Vb(x) and the real and 
imaginary parts of the resonance density of v u (x). For these 
calculations a — 4, b — 0.05, c = 4 and d — 10. Inset (a) 
shows the potential Vb(x) and inset (b) shows the potential 

v„(i) 



imaginary parts correspond respectively to the correct 
energy and inverse-lifetime of the metastable state. 

We note that the lifetime can also be extracted directly 
from the resonance density by looking at its asymptotic 
behavior (r — > oo), just as the ionization energy / of 
bound systems is extracted from the r — > oo decay of 
(real) ground-state densities. For one-electron systems, 
this decay is simply determined by the wavenumber k = 
v2~7- Upon complex scaling, the 'ionization energy' and 
'decay constant' both become complex (J — » Ig, and k — > 



and Ig 



|e Kg. for trie example below, using 
Eq.(j6j, we find Ig = -1/4 + iy/3/4, correctly. 

Next, a numerical example will allow us to look deeper 
into the meaning of the complex density, and how the 
resonance lifetime is extracted from it. 

Numerical Test: We study the fate of the bound state 
of a simple one-dimensional gaussian well, Vb(x) = a(l — 

e~T) ; when steps are imposed on it at id, making the 
problem unbound (see inset of[T]), v u (x): 



y u (x) 
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I _|_ e -2c(x+d) I _|_ g-2c(x-rf) 



. (7) 

We use complex-scaling to solve for the complex eigen- 
values and eigenfunctions. As the steps are taken far 
away [d — > oo), the resonance energy and the real part of 
the resonance density approach the bound state energy 
and ground-state density of the bound potential Vb(x), 
respectively, as expected [251 HH] • 

For example, when a = 4 and b = 0.05 the bound- 
state energy of an electron in Vb{x) is 1.6246. [I] shows 
that this same value is obtained for v u (x) when d > 2. 
Also, the imaginary part of the resonance energy shrinks 
exponentially as the steps get wider. This behavior is 
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TABLE I. Resonance energies for v u (x) with various values 
of d. For all calculations c = 4, a = 4, b — 0.05 . 
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FIG. 2. An example 9 trajectory. In this case the real part of 
the energy is plotted against the imaginary part of the energy 
for 6 values from 0.24 to 0.54, and the optimum 9 is found at 
the kink with a value of 0.312. 



expected because the inverse of QEg corresponds to the 
lifetime of the resonance, and the resonance is long lived 
for wide steps in v u (x). 

The Fourier Grid Hamiltonian technique (FGH) [30, 
|3T| is found to be the most efficient method for these com- 
plex scaling calculations. It allows one to solve the eigen- 
value problem without the use of a basis set expansion. 
The Hamiltonian is written in coordinate representation, 
(x\Hg\x'), and the spatial coordinate is discretized. This 
matrix is then diagonlized and the resulting eigenvectors 
give directly the amplitude of the solutions of the reso- 
nance wave function. Using this method, the resonance 
energy for v^ix) is converged using less than 100 grid 
points. 

The resonance energy and lifetime is independent of 9 
|llj . However when doing a numerical calculation, the 
use of finite basis sets or finite grids necessitates finding 
an optimum 9. This optimum condition is acheived by 
examining ^-trajectories [11] (see [2]). When the real and 
imaginary part of the resonance energy become station- 
ary around the optimum value of 9 one can see a kink 
or loop in the trajectory. Within our range of parame- 
ters, the optimum 9 for v„(x) falls between 0.28 and 0.4. 
With many more grid points this procedure becomes less 
important as all the theta trajectory points collapse to 
the region of the resonance. 

The key point is this: When ng(x) is normalized ac- 
cording to Eq. [2j the correct resonance energies and life- 
times are indeed obtained from e~ 2l0 T\y[ng\. For exam- 
ple, with a = 4, b = 0.05, c = 4 and d = 2 the resonance 
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FIG. 3. The % error in the real and imaginary part of the reso- 
nance energy calculated from Eq|3]with Fe[ne] = e~ 2 ' B Tw[ng] 
as a function of the number of grid point used in the calcu- 
lation. For this plot a — 4, b = 0.05, c = 4 and d = 2. 
Also included is the convergence of the bound energy of v&(a:) 
calculated from the unsealed 2V[nl functional. 



energy is 1.6285-3.3630 x 10 _3 i. Using the FGH method 
with 999 grid points, the von Weizsacker functional gives 
a resonance energy of 1.6282 — 3.1126 x 10 _3 i. The con- 
vergence of the functional calculations is included in [3] 
where the error is defined relative to FGH results ob- 
tained using the same number of grid points. The error 
of the resonance energies obtained with the functional is 
less than 1% for more than 200 grid points, and the error 
of the lifetimes is less than 1% for more than 2500 grid 
points. From [3j it is clear that the real part of the res- 
onance density converges very quickly and its behavior 
is comparable to that of the bound energy convergence 
obtained for Vf,(a;) via The convergence of the 

imaginary part follows a similar trend, but errors of less 
than 10% are only obtained when more than 1000 grid 
points are employed. The imaginary functional is exact, 
but sensitive to small changes in the resonance density. 



CONCLUSION 

For one-electron systems, we have shown that the 
lowest-energy resonance lifetime is encoded into the cor- 
responding complex resonance density, and can be ex- 
tracted from it via the properly scaled density functional. 
If the same result was established for many-electron sys- 
tems, a new way of calculating 'ground-state' resonance 
lifetimes could be based on a complex-scaled version of 
standard KS-DFT, an exciting prospect. We are working 
along these lines. 



ACKNOWLEDGMENTS are gratefully acknowledged. 

Support from the Petroleum Research Fund grant 
No.PRF# 49599-DNI6, and discussions with Yu Zhang, 



[1] R. E. Palmer and P. J. Rous, Rev. Mod. Phys. 64, 383 
(1992) 

[2] J. Simons, J. Phys. Chem. A 112, 6401 (2008) 

[3] R. K. Nesbet, Phys. Rev. A 62, 04070RR) (2000) 

[4] P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964) 

[5] W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 (1965) 

[6] J. C. Rienstra-Kiracofe, G. S. Tschumper, H. F. Schaefer, 

S. Nandi, and G. B. Elliso, Chem. Rev. 102, 231 (2002) 
[7] M. Puiatti, D. M. A. Vera, and A. B. Pierini, Phys. Chem. 

Chem. Phys. 10, 1394 (2008) 
[8] E. Runge and E. K. U. Gross, Phys. Rev. Lett. 52, 997 

(1984) 

[9] A. J. Krueger and N. T. Maitra, Phys. Chem. Chem. 
Phys. 11, 4655 (2009) 
[10] A. Wasserman and N. Moiseyev, Phys. Rev. Lett. 98, 
093003 (2007) 

[11] N. Moiseyev, P. R. Certain, and F. Weinhold, Mol. Phys. 

36, 1613 (1978) 
[12] W. P. Reinhardt, Ann. Rev. Phys. Chem. 33, 223 (1982) 
[13] B. Simon, Ann. Math. 97, 247 (1973) 
[14] E. Balslev and J. M. Combes, Commun. Math. Phys. 22, 

280 (1971) 
[15] Y. K. Ho, Phys. Lett. 99, 1 (1983) 
[16] S. Hein and W. Koch, J. Fluid. Mech. 605, 401 (2008) 
[17] L. I. Deych and I. V. Pnonomarev, Phys. Rev. B 71, 

035342 (2005) 
[18] C. W. McCurdy, Phys. Rev. A 21, 464 (1980) 
[19] B. R. Junker, Adv. At. Mol. Phys. 18, 207 (1982) 
[20] P. Nordlander and J. C. Tully, Phys. Rev. B 42, 5564 

(1990) 

[21] P. Nordlander and J. C. Tully, Surf. Sci. 211, 207 (1989) 
[22] P. Nordlander and J. C. Tully, Phys. Rev. Lett. 61, 990 

(1988) 

[23] D. A. Telnov and S. I. Chu, Phys. Rev. A 66, 043417 
(2002) 

[24] N. Moiseyev, Phys. Lett. 33, 1982 (1982) 

[25] M. Ernzerhof, J. Chem. Phys. 125, 124104 (2006) 

[26] C. F. von Weizsacker, Z. f. Phys. 96, 431 (1935) 

[27] G. D. Doolen, Int. J. Quant. Chem. 14, 523 (1978) 

[28] A. Buchleitner, B. Gremaud, and D. Delande, J. Phys. 

B: At. Mol. Phys. 27, 2663 (1994) 
[29] J. R. Taylor, Scattering Theory (Dover Publications, Inc., 

New York, 1972) 
[30] C. C. Marston and G. G. Balint-Kurti, J. Chem. Phys. 

91, 3571 (1989) 
[31] S. Chu, Chem. Phys. Lett. 167, 155 (1990) 



